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Abstract 

We study the effect of finite chemical potential for the 
QGP constituentsin the Ramanathan et aLstatistical model 
(Phys.Rev.C70, 027903,2004). While the earher computa- 
tions using this model with vanishing chemical potentials 
indicated a weakly first order phase transition for the sys- 
tem in the vicinity of 170 MeV (Pramana, 68, 757, 2007), 
the introduction of finite values for the chemical poten- 
tials of the constituents makes the transition a smooth 
roll over of the phases, while allowing fireball formation 
with radius of a few 'fermi' to take place. This seems to 
be in conformity with the latest consensus on the nature 
of the QGP-Hadron phase transition. 
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1 Introduction 



The formation of QGP droplet (fireball) is one of the most exciting 
possibility in ultra relativistic heavy ion collisions(URHIC) [Ij.The 
physics of such an event is very complicated and to extract meaning- 
ful results from a rigorous use of QCD appropriate to this physical 
system is almost intractable though heroic efforts at lattice estima- 
tion of the problem has been going on for quite some time[2]. One 
way out is to replicate the approximation schemes which have served 
as theoretical tools in understanding equally complicated atomic and 
nuclear systems in atomic and nuclear physics in the context of QGP 
droplet formation. This approach lays no claim to rigour or ab-initio 
understanding of the phenomenon but lays the framework on which 
more rigorous structures may be built depending on the phenomeno- 
logical success of the model as and when testable data emerges from 
ongoing experiments. 

The nucleation process is driven by statistical fluctuations which 
produce the QGP droplets in a hadronic phase,the size of the fluc- 
tuations being determined by the critical free energy difference be- 
tween the two phases. The Kapusta et-al model (3) uses the liquid 
drop model expansion for this, as given by 
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The first term represents the volume contribution,the second 
term is the surface contribution where a is the surface tension, and 
the last term is the so called shape contribution. The shape contribu- 
tion is an entropy term on account of fluctuations in droplet shape 
which we may ignore in the lowest order approximation. The criti- 
cal radius Rc can be obtained by minimising (1) withrespect to the 
droplet radius R, which in the Linde approximation [6] is. 



Rc = — or a 
Ap 



(2) 



2 The statistical model 



In the approximation scheme of the Ramanathan et.al[4 ] the rel- 
ativistic density of states for the quarks and gluons is constructed 
adapting the procedures of the Thomas-Fermi construction of the 
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electronic density of states for complex atoms and the Bethe density 
of states [6] for nucleons in complex nuclei as templates .The expres- 
sion for the density of states for the quarks and gluons(q,g) in this 
model is: 

where k is the relativistic four-momentum of the quarks and the 
gluons,// is the volume of the fire ball taken to be a constant in the 
first approximation and V is a suitable confining potential relevant 
to the current quarks and gluons in the QGP [4] given by: 

V^Qik) = {l/2kh,,, g\k)T-' - ml/2k . (4) 

where m is the mass of the quark which we take as zero for the 
up and down quarks and ISOMev.for the strange quarks. The g(k) is 
the QCD running coupling constant given by 

g\k) = (4/3)(127r/27){l/ln(l + P/A2)} , (5) 

where A is the QCD scale taken to be 150 MeV. 
The model has a natural low energy cut off at: 

kmin = V{k^in) OT k^in = {^i g^^^T^ / 2)^' \ (6) 

with N = [(4/3)(127r/27)]3. 

The free energy of the respective cases i(quarks,gluons,interface 
etc.) for Fermions and Bosons(upper sign or lower sign) can be 
computed using the following expression: 

F, = TTg. J dkp,{k) ln(l ± e'^V^)/^) , (7) 

With the surface free-energy given by a modified Weyl[7] expres- 
sion: 

^^nterfa.e = lTjdkp^,,,ik)6{k-T), (8) 

where the hydrodynamical flow parameter for the surface is: 

7 = X v'(l/7,)2 + (1/^^)2^ (9) 

For the pion which for simplicity represents the hadronic medium 
in which the fire ball resides ,the free energy is: 

= (3r/27r2)zy r k^dk ln(l - g-V-l+^VT) ^ (^g) 
Jo 
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With these ingredients we can compute the free-energy change 
with respect to both the droplet radius and temperature to get a 
physical picture of the fire ball formation, the nucleation rate govern- 
ing the droplet formation ,the nature of the phase transition etc. This 
can be done over a whole range of flow-parameter values [4], We ex- 
hibit only the two most promising scenarios in fig.l and fig. 2. Fur- 
ther investigating the properties of the corresponding free energies 
as in [4], it was found that the QGP-hadron phase transition is a 
weakly first order one. 
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Figure 1: Ffotai at 7^ = 67^, 7g = 1/6 for various temperatures. 
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Figure 2: Ftotai at Ig = 87^ , 7g = 1/6 for various temperatures. 
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3 Effect of finite chemical potential on the phase 
transition 

All our previous calculations were based on the assumption of van- 
ishing chemical potential for the quarks, which is not a realistic one. 
Therefore, we consider the changes in our result by introducing fi- 
nite chemical potentials into our formulation. This is achieved by ,as 
is done in statistical mechanics ,by changing the exponents in (7), 
from — ^ to where '//'is the chemical potential. Following the 
typical values of the chemical potentials for the quarks and gluons in 
the literature [8], we use at 300 MeV and 400 MeV for the different 
parameter values of the model. The result of our computations are 
displayed in figs. 3 to 14. 
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Figure 3: Variation of S with temperature T at 7^ = 67^ , 7^ = 1/6. 
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Figure 4: Variation of S with temperature T at 7^ = 87, , 7^ = 1/6. 
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Figure 5: Variation of specific fieat Cy witli temperature T at 7g = 

679 , 7g = 1/6- 
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Figure 6: Variation of specific fieat Cy witli temperature T at 7g = 

, Iq = 1/6- 
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Figure 8: Variation of S with temperature T at 7^, = 87^ , 7^ 
at /i = 300 . 
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Figure 9: Variation of S with temperature T at 7^, = 67^ , 7^ = 1/6 
at /i = 400 . 



14 




150. 160 170 18,0. 190 



T(Me\/) ATr=8rand/i^00 



Figure 10: Variation of S with temperature T at 7^ = 87^ , 7^ = 1/6 
at /i = 300 . 
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Figure 11: Variation of Cy with temperature T at 7^ = 67^ , % = 1/6 
at /i = 300 . 
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Figure 14: Variation of Cy with temperature T at 7^ = 87^ , % = 1/6 
at = 400 . 
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As could be seen from the above graphs the effect of finite chem- 
ical potential seems to make the phase transition a smooth roll over 
from the hadron phase to the QGP phase rather than a weakly first 
order transition predicted by the model for the unrealistic zero chem- 
ical potential case [4]. This result is in consonance with the present 
expectations from lattice simulations [8]. 
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